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ABSTRACT 

We present a dynamical model of the high mass X-ray binary LMC X-1 
based on high-resolution optical spectroscopy and extensive optical and near- 
infrared photometry. From our new optical data we find an orbital period of 
P = 3.90917±0.00005 days. We present a refined analysis of the All Sky Monitor 
data from RXTE and find an X-ray period of P = 3.9094 ±0.0008 days, which is 
consistent with the optical period. A simple model of Thomson scattering in the 
stellar wind can account for the modulation seen in the X-ray hght curves. The 
V - K color of the star (1.17 ± 0.05) implies Ay = 2.28 ± 0.06, which is much 
larger than previously assumed. For the secondary star, we measure a radius of 
i?2 = 17.0±0.8 Rq and a projected rotational velocity of V^rot sin i — 129.9±2.2 km 
s~^. Using these measured properties to constrain the dynamical model, we find 
an inclination ofi = 36.38±1.92°, a secondary star mass of M2 = 31.79±3.48 Mq, 
and a black hole mass of 10.91 ± 1.41 Mq. The present location of the secondary 
star in a temperature-luminosity diagram is consistent with that of a star with 
an initial mass of 35 Mq that is 5 Myr past the zero-age main sequence. The 
star nearly fills its Roche lobe (~ 90% or more), and owing to the rapid change 
in radius with time in its present evolutionary state, it will encounter its Roche 
lobe and begin rapid and possibly unstable mass transfer on a timescale of a few 
hundred thousand years. 



1. Introduction 

The first X-ray source to be discovered in the Magellanic Clouds, LMC X-1 (Mark et al. 
1969), is a persistently luminous (L^ > 10^^ erg s~^ ) X-ray binary that has been observed 
by nearly all X-ray missions during the past 37 years (e.g., Leong et al. 1971; Cui et al. 
2002). Spectroscopic studies of its optical counterpart revealed an orbital period of 4 days 
and a probable mass for the compact star "near M ^ 6 Mq" (Hutchings et al. 1983, 1987), 
making it the fourth dynamical black-hole candidate to be established. 

The candidacy of this black hole has had a checkered history. Even the most fundamental 
property of this binary system - its orbital period - was established only very recently. For 

nearly 20 years, the accepted orbital period was 4.2288 d, the value adopted in the dynamical 
study of Hutchings et al. (1987). The correct orbital period, P = 3.9081 ± 0.0015 d, was 
determined by Levine & Corbet (2006) using i^XTE" All-Sky Monitor (ASM) X-ray data. The 
optical results reported in this paper and further analysis of the ASM data amply confirm 
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that the Levine & Corbet (2006) period is correct. As an interesting side note, the correct 
period is noted in the earher Hutchings et al. (1983) work as one of several candidate orbital 
periods. 

Even the identity of the optical counterpart was not firmly established until fairly re- 
cently. Initially, based on a rather uncertain X-ray position, a B5 supergiant known as R148 
was favored over what is now known to be the counterpart, an 07/08 giant identified as 
"star 32" by Cowley, Crampton & Hutchings (1978). The counterpart was finally established 
through an analysis of multiple ROSAT HRI observations by Cowley et al. (1995). A precise 
Chandra X-ray position (Cui et al. 2002) and the agreement between the X-ray and optical 
periods (mentioned above) leave no doubt that star 32 of Cowley et al. is the counterpart of 
LMC X-1. We present high resolution V- and i^-band finding charts of the field in Figured! 

In this paper, we confirm the basic model and the principal conclusions presented by 
Hutchings et al. (1987) while greatly improving upon their pioneering work. Thirty high 
resolution and 14 medium-resolution spectra have allowed us to reduce the uncertainty in 
the radial velocity amplitude by a factor of 6 and to obtain a secure value for the rotational 
line broadening. Of equal importance, we present the first optical light curves and infrared 
magnitudes and colors of LMC X-1, the analysis of which allows us to strongly constrain 
the orbital inclination angle and hence our dynamical model of the binary. In ^ and ^ 
we present our new photometric and spectroscopic observations. In ^we present improved 
measurements of the orbital period, the temperature of the secondary star, and the radius 
of the secondary star. These new observations are used to construct a dynamical model of 
the system, which is presented in ^ a summary of our results is in ^ 

2. Photometry 

2.1. Observations and Reductions 

Optical and near infrared observations of LMC X-1 were obtained on 66 nights between 
2007 January 6 and 2007 April 2 using the ANDICA1V0 instrument on the 1.3m telescope at 
CTIO, which is operated by the Small and Moderate Aperture Research Telescope System 
(SMARTS) consortiunll. During each visit to the source, a 120 second exposure in the 
B filter and a 120 second exposure in the V filter were obtained simultaneously with an 
exposure in J that consisted of 10 dithered subexposures of 12 seconds each. Flat field 
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images for all filters were taken nightly, and IRAF tasks were used to perform the standard 
image reductions (i.e. bias subtraction and flat-fielding for the optical and fiat-fielding, sky 
subtraction, and image combination for the infrared). All of the images were inspected 
visually, and a few of them were discarded owing to poor signal-to-noise. We retained 58 
images in B, 57 images in V, and 57 images in J. 

We also observed the field of LMC X-1 on the nights of 2006 December 6 and December 
8 using the 6.5m Magellan Baade telescope at Las Campanas Observatory (LCO) and the 
Persson's Auxiliary Nasmyth Infrared Camera (PANIC; Martini et al. 2004). On December 
6 three composite frames were obtained in the J band and four composite frames were 
obtained in the K band in 0'.'7 seeing. Each composite frame consists of five 7-second 
dithered subexposures. On December 8, ten sequences in J with 3-second subexposures and 
eight sequences in K with 7-second subexposures were obtained in 0'.'5 seeing. The conditions 
were photometric on both nights. The PANIC data were reduced and processed with IRAF 
and custom PANIC software to produce mosaic frames of the field (13 mosaic frames in J 
and 12 mosaic frames in K). 



2.2. Derivation of the Photometric Light Curves 

PSF-fitting photometry was used on the PANIC mosaic images to obtain instrumental 
magnitudes of LMC X-1 and two nearby stars. The absolute calibration was done with 
respect to stars from the 2MASS catalog (Skrutskie et al. 2006). Comparison stars having 
less than 8000 detector counts in the PANIC frames were selected and checked for variability. 
Weighted differential photometry of LMC X-1 was performed with respect to seven and 
six comparison stars for the J and K-hand mosaic images, respectively. For LMC X-1, 
we found a mean J magnitude of J = 13.76 with an rms of 0.02 mag and a mean K 
magnitude of K = 13.43 with an rms of 0.01 mag. Optimal aperture photometry was also 
performed for LMC X-1 and the comparison stars, and we obtained similar values for the 
mean magnitudes. For comparison, the magnitudes for LMC X-1 given in the 2MASS catalog 
are J = 13.695 ± 0.063 and K = 13.293 ± 0.063, although these values should be treated 
with caution owing to the presence of the bright star R148 ~ 6" to the east (Figured]). 

The photometric time series for the B and V filters were obtained using the programs 
DAOPHOT He, ALLSTAR, and DAOMASTER (Stetson 1987, 1992a, 1992b; Stetson, Davis, 
& Crabtree 1991). The instrumental magnitudes were placed on the standard scales using 
observations of the Landolt (1992) fields RU 149, PC 1047+003, and PC 1657+078. Aperture 
photometry was done on the standard stars and on several bright and isolated stars in the 
LMC X-1 field and the DAOGROW algorithm (Stetson 1990) was used to obtain optimal 
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magnitudes. We find for LMC X-1 an average V magnitude of ^ = 14.60±0.02, and a, B — V 
color of i? — = 0.17 ± 0.08, marginally consistent with the values of ^ = 14.52 ± 0.05 and 
B -V = 0.29 ± 0.02 given in Bianchi & PakuU (1985). 

The photometric time series for the SMARTS J band was obtained using aperture 
photometry. A relatively small aperture radius (~ 1") was used to exclude light from nearby 
sources. Stars from the 2MASS catalog were used to place the instrumental magnitudes 
onto the standard system. The J-band magnitudes from SMARTS are in agreement with 
the mean J-band magnitude derived from the PANIC data. 

Figure [2] shows the light curves phased on the photometric ephemeris determined below. 
The light curves show the double-wave modulation characteristic of ellipsoidal variations, 
with maxima at the quadrature phases and minima at the conjunction phases. The amplitude 
of the modulation [k. 0.06 mag, maximum to minimum) is not especially large, which usually 
indicates a relatively low inclination, or a small Roche lobe filling factor for the mass donor, 
or both. 



3. Spectroscopy 

3.1. Observations and Reductions 

Thirty spectra of LMC X-1 were obtained on the nights of 2005 January 19-24 using 
the Magellan Inamori Kyocera Echelle (MIKE) spectrograph (Bernstein et al. 2002) and 
the 6.5 m Magellan Clay telescope at LCO. The instrument was used in the standard dual- 
beam mode with a I'.'O x 5'.'0 slit. Exposure times ranged between 1200 and 2700 seconds 
and observing conditions were good. The seeing was well below 1" most of the time with 
excursions as low as 0'.'5 and as high as 1'.'2. The pair of 2048 x 4096 pixel CCD detectors 
were operated in the 2x2 on-chip binning mode. The blue arm had a wavelength coverage 
of 3340-5065 A and the spectral dispersion on the MIT detector was A/AA = 100, 000 (0.03- 
0.05 A pixel"^), whereas the red arm covered 4855-9420 A and the dispersion on its SITe 
detector was A/AA = 71, 000 (0.07-0.13 A pixel^^). Our 1" slit delivered a spectral resolution 
of R = 33, 000 and 28,000 in the blue and red arms respectively. ThAr lamp exposures were 
obtained before and after each pair of observations of the object, and several fiux standards 
and spectral-comparison stars were observed. 

MIKE was located at a Nasmyth focus, which minimizes fiexure and calibration prob- 
lems. No instrument rotator or dispersion compensation optics were available, and the 
position angle of the spectrograph was set to the parallactic angle for an object at air mass 
1.3. Because all of our observations were taken below air mass ~ 1.5, the effects of light loss 
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due to atmospheric dispersion are small and unlikely to significantly affect our results. 

The MIKE data were reduced using a pipeline written by Dan Kelsor§|. The pipeline 
performs all of the detector calibrations in the standard way using bias frames and tungsten 
fiat fields which were obtained each afternoon. The arc spectra talcen at the position of each 
target were used to correct for the non-orthogonality between the dispersion and spatial axes 
and construct a 2D wavelength solution for all object frames. All object orders were then 
sky subtracted using the technique discussed in Kelson (2003) and optimally extracted. The 
signal-to-noise ratio per pixel in the orders with useful lines was generally in the range of 
~ 50 — 100 for most of the spectra. 

The extracted spectra were inspected visually and artifacts due to cosmic rays were 
removed manually using simple interpolation. In addition to cosmic rays, there were some 
artifacts in the cores of the higher Balmer lines owing to imperfect subtraction of nebular 
lines. These artifacts were also removed by interpolation. 

Finally, the individual orders in the cleaned spectra were trimmed to remove the low 
signal-to-noise regions at each end and normalized using cubic spline fits. The trimmed and 
normalized orders were then merged into a pair of spectra (e.g. red arm and blue arm). 
Figure [3] shows the average spectrum (in the restframe of the secondary) for much of the 
blue arm and part of the red arm. The model spectrum shown in Figure [3] is discussed in 

An additional fourteen spectra of LMC X-1 were obtained on the nights of 2008 August 
2-3 using the Magellan Echellette Spectrograph (MagE; Marshall et al. 2008) and the Clay 
telescope at LCO. We used the 0'.'7 slit which yields a resolving power oi R = 6000. The 
exposure times were 600 seconds and the seeing ranged from 0'.'8 to I'.'l for the first night 
and 1'.'4 to 2'! 7 for the second night. ThAr lamp exposures were obtained before and after 
each pair of observations of LMC X-1. 

The images were reduced with tasks in the IRAF 'ccdproc' and 'echelle' packages. After 
the bias was subtracted from each image, pairs of images were combined using a clipping 
algorithm to remove cosmic rays. The resulting seven images were fiat-fielded using a nor- 
malized master fiat and then rotated in order to align the background night sky emission 
lines along the columns. After the background emission lines were rectified, the spectra from 
individual orders were optimally subtracted. Unfortunately, the rectification of the back- 
ground lines was not perfect, which produced artifacts due to nebular lines in the cores of 
the higher Balmer lines. These imperfections were removed by simple interpolation. The 
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process used to merge the orders in the MIKE spectra described above was also used to 
merge the MagE orders. The resulting seven spectra cover a useful wavelength range of 
3200-9250 A. 



3.2. Radial Velocities 

In a recent study of three high mass X-ray binaries with OB-star companions, van der 
Meer et al. (2007) measured radial velocity curves using single hues. They found that X-ray 
heating produced distortions in some of the line profiles, thereby yielding for some cases 
different X-velocities for different lines. 

In order to assess the effects of X-ray heating on our spectra, we selected four bandpasses 
with strong lines that have a wide range of excitation energies over which to derive radial 

velocities from the MIKE spectra: 3780-3880 A, dominated by H9 & HIO and including He 
I A3820; 3990-4062 A, He I A4026; 4150-4250 A, He II A4200; and 4425-4525 A, He I A 4471. 
These lines were selected because of their prominence in both object and template spectra 
and their isolation from other spectral features. 

The radial velocities were determined using the FXCOR task within IRAF, which is an 
implementation of the cross-correlation analysis developed by Tonry & Davis (1979). For 
each bandpass, the spectra were low-frequency filtered either by fitting a Legendre polynomial 
or by Fourier filtering. To optimize the removal of any order blaze remnants, the Legendre 
polynomial orders and the Fourier ramp filter components were determined individually 
for each bandpass. For example, the region containing H9 and HIO required a Legendere 
polynomial order of 8 (i.e. nine terms) or a Fourier cut-on wavenumber of 4 and a full- 
on wavenumber of 8, whereas the region containing the fine He I A4471 needed a 28-piece 
polynomial or a cut-on of 9 and a full-on of 18. The spectra were then continuum-subtracted 
prior to computing the cross-correlation. No other filtering was applied. 

The correlations were performed using three different templates: a model spectrum from 
the OSTAR2002 grid (Lanz & Hubeny 2003; see the discussion below) and the spectra of 
two of our comparison stars: HD 93843, 06 Ill(f) and HD 101205, 07 IIIn((f)) (Walborn 
1972, 1973). Typical values of the Tonry & Davis (1979) r-value, a measure of signal-to-noise 
ratio, ranged from about 10 to 70 and the median value was 26. 

The radial velocity data for each bandpass were fitted to a circular orbit model which 
returns the systemic velocity Vq, the time of maximum velocity Tmax, and the velocity semi- 
amplitude of the secondary K2- The orbital period was fixed at our adopted value, which 
is given in Table 3. The model provided good fits to the data. For each fit, the statistical 
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errors on the velocities returned by the fit were scaled by factors ranging from 0.3 to 2.1 and 
the data were refitted in order to give a reduced chi-squared {xD of unity. 

As an aid to visualize the systematic effects, Figure H] summarizes the fitted values of K2 
and Tinax obtained for different combinations of (i) the four bandpasses, (ii) the three template 
stars and (iii) the two modes of low-frequency filtering, Legendre and Fourier. The open 
circles show, for example, the results obtained by fitting the velocity data using the spectrum 
of HD 93843 as the template and employing Legendre filtering. The other open symbols show 
for comparison the results obtained with the other template/filter combinations. The results 
shown as filled circles and bold error bars were derived using the weighted mean velocities 
for each spectrum: K2 = 69.19±0.88 km s"^ and T^ax = HJD 2453392.3043±0.0081. These 
mean values are indicated by the dashed lines. 

Figure S] shows that our results for K and T^ax are quite insensitive to the choice of 
bandpass, template or mode of filtering, although there are some systematic effects at play. 
For example, the weighted mean K-velocitj of the He II line, ii'(A4200)= 70.87 ± 0.81 km 
s~^, is greater than our adopted mean value, albeit still within 2a. This modest shift in K is 
far less than was observed in the recent van der Meer et al. (2007) study of three 0-star X-ray 
binaries that contain neutron star primaries. Because we have considered lines with a wide 
range of excitation energies (Balmer, He I and He II) and have found quite consistent results, 
we conclude that our dynamical results are little affected by X-ray heating, tidal distortion 
of the secondary, stellar wind, and other commonly-observed sources of systematic effects. 

Finally, we measured the equivalent widths of the HIO, He I A4471, and He II A4200 
lines in the individual spectra. There is no apparent trend with orbital phase (see Figure [5]), 
which further confirms that the effects of X-ray heating are minimal. 

Knowing that the velocities are little effected by X-ray heating, we derived the final 
adopted velocities from the MIKE and MagE spectra using FXCOR, a synthetic template, 
and a cross correlation region that covers nine He I and He II lines between 4000 A and 
5020 A. The mean velocity curve is shown plotted in Figure [61 

4. Results 

4.1. Refined Orbital Period 

Levine & Corbet (2006) analyzed RXTE ASM data of LMC X-1 taken between 1996 
March and 2005 November. They reported a periodicity of 3.9081 ± 0.0015 days, which 
differs from the period of P = 4.2288 days given in Hutchings et al. (1987). On the other 
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hand, the X-ray period is consistent with one of the possible periods given in Hutchings et 
al. (1983), namely P = 3.909 ± 0.001 days. 

We derived a refined period from our optical data. The optical light curves completely 
rule out the 4.2288 day period and strongly favor the period near 3.909 days. We made a 
periodogram from the Magellan radial velocities by fitting a three-parameter sinusoid to the 
data at various trial periods and recording the of the fit. The results are shown in the top 
of Figure [71 A unique period cannot be found from these data alone, leaving us with possible 
periods near 3.884, 3.898 and 3.909 days. Although the radial velocities given in Hutchings 
et al. (1983, 1987) have relatively large errors and the times given are accurate to only about 
one minute, they can be used to rule out possible alias periods. When these additional 
velocities are included in the analysis, only the period near 3.909 days remains viable (see 
the center panel in Figure [T]). Finally, we used the ELC code (Orosz & Hauschildt 2000) to 
compute periodograms by modeling both our light and velocity curves together. The details 
of the model fitting are discussed thoroughly in §4.4.11 For the purposes of computing 
a periodogram, the period is held fixed at a given trial period and several other fitting 
parameters are varied (here the orbit is assumed to be circular) until the fit is optimized. 
The graph of vs. P (see the bottom panel of Figure [7]) then serves as a periodogram. The 
combination of the SMARTS photometry and the Magellan radial velocities yields a unique 
period of 3.90917 ± 0.00005 days. This period is consistent with the period found from the 
combination of the Hutchings et al. (1983, 1987) and Magellan radial velocities. This period 
is also consistent with the X-ray period found by Levine & Corbet (2006), and with the 
refined X-ray period that we now derive. 

The ASM c onsists of three Sca nning Shadow Cameras (SSCs) mounted on a rotating 
Drive Assembly ( Levine et al.lE996l ). Approximately 53,700 measurements of the intensity 



of LMC X-1, each from a 90-s exposure with a single SSC, were obtained from the beginning 
of the RXTE mission in early 1996 through 2008 June. A single exposure yields intensity 
estimates in each of three spectral bands which nominally correspond to photon energy 
ranges of 1.5-3, 3-5, and 5-12 keV with a sensitivity of a few SSC counts s~^ (the Crab 
Nebula produces intensities of 27, 23, and 25 SSC counts s~^ in the 3 bands, respectively). 
The X-ray intensity of LMC X-1, as seen over more than 12 years in the ASM light curves, 
has been more or less steady near 20 mCrab (1.5-12 keV) with variations, when 10-day 
averages are considered, of ~ ±10%. 

A periodicity in the X-ray intensity of LMC X-1 was found during a search for period- 
icities in ASM data using advanced analysis techniques. Two strategies were applied in this 
search to improve the sensitivity. The first involves the use of appropriate weights such as 
the reciprocals of the variances, since the individual ASM measurements have a wide range 
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of associated uncertainties. The second strategy stems from the fact that the observations 
of the source are obtained with a low duty cycle, i.e., the window function is sparse (and 
complex). The properties of the window function, in combination with the presence of slow 
variations of the source intensity act to hinder the detection of variations on short time 
scales. The power density spectrum (PDS) of the window function has substantial power 
at high frequencies, e.g., 1 cycle d~^ and 1 cycle per spacecraft orbit (~ 95 minute period). 
Since the data may be regarded as the product of a (hypothetical) continuous set of source 
intensity measurements with the window function, a Fourier transform of the data is equiva- 
lent to the convolution of a transform of a continuous set of intensity measurements with the 
window function transform. The high frequency structure in the window function transform 
acts to spread power at low frequencies in the source intensity to high frequencies in the 
calculated transform (or, equivalently, the PDS). This effectively raises the noise level at 
high frequencies. 

In our analysis, the sensitivity to high frequency variations is enhanced by subtracting 
a smoothed version of the light curve from the unsmoothed light curve. To perform the 
smoothing, we do not simply convolve a box function with the binned light curve data since 
that would not yield any improvement in the noise level at high frequencies. Rather, we 
ignore bins which do not contain any actual measurements and we use weights based on 
estimates of the uncertainties in the individual measurements to compute the smoothed 
light curve. The kernel used in the smoothing was a Gaussian with a full width at half 
maximum of 30.0 days, so the smoothed light curve displayed only that variability with 
Fourier components at frequencies below ^ 0.03 day^^. The smoothed light curve was 
subtracted from the unsmoothed light curve, and the difference light curve was Fourier 
transformed. The results are illustrated in Figure [81 The center frequency of the peak 
and the formal uncertainty thereof are 0.25580 ± 0.00005 d~^. The corresponding period is 
3.9094 ± 0.0008. This period differs from the optical period by only 1.3 x 10~^ days. 

In Figure [9] we show folded X-ray light curves for the 1.5-3, 3-5, and 5-12 keV photon- 
energy bands as well as for the overall 1.5-12 keV band. There is a clear orbital modulation 
with the minimum intensity occurring at phase 0, which corresponds to the time of the 
inferior conjunction of the secondary star. The X-ray light curves can be fitted to a function 
of the form /(0) = + ai cos(27r0) -|- a2 sin(27r0). Table [T] gives the best-fitting coefficients 
for the three individual bands and for the sum. In all cases, the a2 terms are consistent with 
zero. We discuss this modulation in more detail in §5.31 
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4.2. Temperature and Rotational Velocity 

The effective temperature of the secondary star and its (projected) rotational velocity 
are needed as inputs and constraints on the dynamical model discussed below. Previously, 
Hutchings et al. (1983) derived a spectral type of 07 for the secondary star based on the 
line strengths of He I, He II, Si IV, and Mg II and a rotational velocity of VT-otsini ^ 150 
km s~^. Negueruela & Coe (2002) derived a spectral type of 08III based on the ratio of the 
He II A4541 and He I A4471 lines and on the strength of the Mg II A4481 line. A spectral 
type of 08III corresponds to an effective temperature of around 33,000 K (Heap, Lanz, & 
Hubeny 2006). 

We used model spectra from the OSTAR2002 grid (Lanz & Hubeny 2003) to derive im- 
proved values of the temperature and rotational velocity of the secondary. The temperature 
and gravity sampling of the grid was extended to 50 K and 0.05 dex, respectively, by first 
resampling the high dispersion model spectra to a fine wavelength spacing common to all 
models and then by interpolation of the resampled models. Hubeny's program R0TINS3 
was used to broaden model spectra with instrumental and rotational broadening profiles, and 
to resample the resulting models to match the wavelength sampling of the observed spectra. 

The "surface" gravity of the secondary star is tightly constrained by the dynamical 
solution discussed below and is found to be logg = 3.50 with a formal la error less than 0.02 
dex. For the purposes of finding the temperature we allow logg to fall within the generous 
range of 3.4 < logg < 3.6. With our large grid of model spectra over a large range of 
temperatures, gravities, and rotational velocities, we performed tests using nine helium 
lines between 4000 A and 5020 A. The line ratios of He II A4541/He I A4471 and He II 
A4200/He I(+II) A4026, which are the main classification criteria for 0-stars (e.g. Walborn 
& Fitzpatrick 1990), are included in this range. We find a temperature of T2 = 33, 225 ±75 K 
from the average MIKE spectrum and T2 = 33, 200±150 K from the average MagE spectrum. 
Figure [3] shows the best-fitting model with the MIKE spectrum. Generally the lines in 
the model spectrum match up with the lines in the observed spectrum. There are a few 
exceptions such as the O II + C HI blend near 4035 A, the Bowen blend near 4640 A, the He 
II line near 4686 A, and the C HI line near 5690 A. The wings and the cores of the Balmer 
lines are not well fit. However, in the case of the MIKE spectra it is difficult to normalize the 
profiles of the Balmer lines owing to their large widths which can take up a substantial part 
of an individual order. The Balmer lines are matched better in the MagE spectrum where 
it is much easier to normalize the individual orders owing to the lower spectral resolution. 

The errors on the temperature quoted above are the formal la errors. The true un- 
certainty on the effective temperature is no doubt larger owing to the fact that the model 
spectrum is not a perfect fit to the data. Also, there are almost certainly systematic errors in 
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the model atmospheres as well. Accordingly, we will adopt a temperature of T2 = 33, 200 K 
with an uncertainty of 500 K for the analysis presented below. 

To refine our initial measurement of the rotational velocity of the star, we used the 
He I lines near 4026.0 A and 4471.0 A and the He II lines near 4541.4 A and 5411.4 A. We 
computed large grids of models with various temperatures, gravities, and rotational velocities 
with the appropriate resolving power for the MIKE spectrograph and performed tests to 
the composite line profiles (here we assume that the main broadening mechanism in the line 
profiles is rotation). Figure fTOl shows the results. The model profiles fit the observed profiles 
fairly well. The rotational velocities found from the individual lines range from 127.9 ± 1.0 
km s~^ to 133.2 ±1.2 km s~^, where the uncertainties are the formal la errors. The average 
value is 129.9 km s~^, and the standard deviation of the four measurements is 2.18 km s^^. 
Hence for the analysis below, we adopt a value of V^ot sini = 129.9 ± 2.2 km s~^. 

4.3. Radius of the Secondary Star and Reddening 

Unlike the case for most Galactic X-ray binaries, the distance to LMC X-1 is well- 
determined. Given the distance, we can find the radius of the secondary star if we can 
obtain good measurements of the apparent magnitude of the star, its effective temperature, 
the extinction, and the bolometric correction. If the radius of the secondary star can be 
found independently, then the dynamical model of the system is further constrained. 

The apparent magnitude of the system is easy to measure. The spectral type of the 
star yields its effective temperature, and the bolometric corrections can be estimated using 
detailed model atmosphere computations (e.g. Lanz & Hubeny 2003; Martins & Plez 2006). 
The remaining quantity, the extinction to the source, is unfortunately subject to the largest 
uncertainties and difficult to measure. 

Using spectra from the International Ultraviolet Explorer (lUE), Bianchi & Pakull (1985) 
determined a color excess of E{B — V) = 0.37 (0.32 from extinction in the LMC and 0.05 
Galactic foreground extinction). Bianchi & Pakull (1985) noted that the color excess derived 
from the lUE data is lower than the value derived from the optical color of the star. Their 
value ofB-V = 0.29 and {B - V)o = -0.27 (Martins & Plez 2006) gives E{B -V) = 0.56 
Using our measured value of S — V = 0.17 we find E{B — V) = 0.44. Thus we have a range 
of color excess values of 0.37 < E{B -V) < 0.56. Using Ay = RvE{B - V) with Ry = 3.1 
(the mean value for the LMC, Cardelli, Clayton, & Mathis 1989), the V^-band extinction 
would be in the range of 1.15 < Ay < 1.74. 

These values of Ay < 1.74 mag are at odds with the extinction derived from two other 
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techniques. In the first case, from our photometry we find a V — K color of V — K = 
1.17 ± 0.05. 0-stars, including dwarfs, giants, and supergiants of 0-subtypes 6 through 9, 
have V — K colors in a relatively narrow range of —0.87 < V — K < —0.83 (Martins & Plez 
2006). Assuming Ak = O.llAy (Cardelh, Clayton, & Mathis 1989), we need Ay = 2.28±0.06 
to match the observed V — K color. In the second case, the \^-band extinction can be 
inferred from the hydrogen column densities derived from X-ray spectra. Table [2] gives seven 
determinations of the column density Nh and the inferred V^-band extinction, assuming 
Av = A^iy/(1.79±0.03) X lO^i, where Nh is given in units of cm'^ (Predehl & Schmitt 1995). 
The values of the extinction cover a considerable range: 2.57 ±0.12 < Ay < 4.53 ± 0.29. 
We believe the first measurement given in the table {Ay = 2.57 ± 0.12 from Cui et al. 2002) 
is the most reliable given the model used (the mult-temperature disk blackbody model) and 
the data (high resolution grating spectra from Chandra HETG). This measurement differs 
by the extinction derived from the V — K color by only 2.5a. Furthermore, we note that the 
X-ray column density may systematically exceed the value inferred from optical reddening 
(i.e. dust) considerations, if even a fraction of the dense stellar wind from the 0-star is not 
completely ionized. 

Adopting Ay = 2.28 mag, the range of color excess values quoted above (0.37 < E{B — 
V) < 0.56) gives 4.07 < Ry < 6.16. Our measured value of color excess of E{B — V) = 0.44 
gives Ry = 5.18. The nominal average value in the LMC is Ry = 3.1, although certain lines- 
of-sight can have Ry values that differ significantly from the average (Cardelli, Clayton, & 
Mathis 1989). 

To compute the radius of the star, we assume a distance modulus to the LMC of 
18.41 ± 0.10 (see Orosz et al. 2007). Although the uncertainty of the distance modulus 
to the center of the LMC itself is somewhat smaller than 0.1 mag, we adopt an uncertainty 
of 0.1 mag to account for the unknown relative position of LMC X-1 with respect to the LMC 
center. We use bolometric corrections for the V and K bands derived from the OSTAR2002 
models with the LMC mean metalicity (Lanz & Hubeny 2003)0 We use apparent magnitudes 
of V = 14.6 ± 0.05 and K = 13.43 ± 0.03 (where we have increased the uncertainties to be 
conservative), and assume the star is the only source of optical light in the system (see the 
discussion in the Appendix regarding the likely faintness of the accretion disk). We use an 
effective temperature of T2 = 33, 200 ± 500 K, and set the gravity to log (7 = 3.50 ± 0.02, as 
determined from the dynamical model described below. We computed the stellar radius using 
the V"-band magnitudes and K-band magnitudes using a wide range of values of Ay, and 



''The OSTAR2002 bolometric corrections are in good agreement with the bolometric corrections given in 
Martins & Plez (2006). We prefer the former since they are tabulated in terms of the effective temperature 
and gravity, although the bolometric corrections are insensitive to the gravity. 
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the results are shown in Figure [TT] The hues cross near Ay ^ 2.2 mag (where R2 ~ 17 Rq), 
which is not at all surprising since we derived Ay = 2.28 ± 0.06 mag from the V — K color 
above. Note that the radius derived from the K band is relatively insensitive to the assumed 
value of Ay, owing to the fact that A^ = 0.11 Ay. The radii found using the K-hand are 
in the range 15.5 ^ R2 ^ 18.3 Rq (this includes the la errors) for extinction values in the 
generous range of 1.4 < Ay < 2.8 mag. For a given value of Ay, the uncertainty in the 
radius derived from the V band is on the order of 1.0 Rq when the generous uncertainties 
in the apparent magnitudes and the large range in effective temperatures are used. For the 
K band, the corresponding uncertainty is on the order of 0.8 Rq. For the purposes of the 
dynamical model below we adopt R2 = 17.0 ± 0.8 Rq, which is the value we find using the 
i^'-band and Ay = 2.28 ± 0.06. The luminosity of the star is then log L/Lq = 5.50 ± 0.05. 

As noted above, the extinction derived from the lUE spectra is unusually low. In order 
to see if the lUE spectra could be explained by a different extinction law, we obtained the 
lUE spectra of LMC X-1 from the MAST archive. The spectra have been reprocessed using 
the latest calibrations, and we combined all of the observations to increase the signal-to- 
noise. The mean lUE spectrum was dereddened using Ay = 2.28 and Ry = 5.18. The 
results are shown in Figure [T^l We also show a model spectrum from the OSTAR2002 grid 
with Teff = 32, 500 K and log (7 = 3.5, scaled using a distance modulus of 18.41 and a stellar 
radius of 17.0 i?©. The agreement between the data and the model is good redward of the 
2200 A bump. However, the 2200 A bump itself seems over-corrected, and the slope of 
the dereddened spectrum blueward of the bump is much too flat. However, the following 
caveats should be noted regarding the lUE data. First, the signal-to-noise of the spectra is 
not terribly high (on the order of 10 per pixel or less). Second, LMC X-1 was observed in 
the lUE large aperture (10"), and the source was not centered in the aperture in order to 
exclude light from the much brighter star R148. Finally, PakuU & Angebault (1986) showed 
that LMC X-1 is in a compact He HI nebula, which is in turn embedded in a larger H II 
region known as N159 (Bianchi & Pakull 1985). Considering the very large aperture of the 
lUE spectra, the UV spectrum of LMC X-1 might be contaminated by scattered light from 
nearby bright stars. It would be worthwhile to obtain additional UV spectra of LMC X-1 to 
further investigate the UV extinction. 

4.4. System Parameters 

4.4- i- Simultaneous Fits to Light and Velocity Curves 

We have several "observables" of the LMC X-1 binary system (e.g. the radial velocity 
curve of the secondary, its B and V light curves, etc.), and we seek the physical model 
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whose observed properties best match the observed data. While the "forward" problem 
of computing the observable properties of a binary system is relatively straightforward, the 
"inverse" problem of deriving a physical model from observed data is much more challenging. 
The ELC code (Orosz & Hauschildt 2000) is a comprehensive code for computing the forward 
problem using a model based on Roche geometry. ELC can also solve the inverse problem 
using its associated optimizing codes based on various numerical techniques (e.g. a "grid 
search", an "amoeba", a Levenberg-Marquardt scheme, a genetic code, and a Monte Carlo 
Markov Chain scheme). 

The ELC model as applied to LMC X-1 has several parameters, many of which can be 
set to reasonable values based on various observed properties of the system. The important 
parameters are related to either the geometrical properties of the system or to the radiative 
properties of the system. 

The free parameters which control the basic system geometry include the orbital period 
P, the orbital separation a (which then gives the total mass via Kepler's third law), the 
ratio of the masses Q = M/M2, the inclination i, and the Roche lobe filling factor /2 of the 
secondary star. We find that the search of parameter space is made more efficient if the 
i^'-velocity of the secondary star K2 and the mass of the secondary star M2 are used to set 
the system scale instead of Q and a. When modeling actual data, one must also specify the 
phase zero point Tq, which is either the time of the inferior conjunction of the secondary star 
(in the case of a circular orbit) or the time of periastron passage (in the case of an eccentric 
orbit). We initially assume the secondary star is rotating synchronously with its orbit, and 
that the star's rotational axis is perpendicular to the orbital plane. We discuss in §4.4.2l how 
our results depend on the assumption of synchronous rotation. 

The parameters which control the radiative properties include the average temperature 
of the secondary star T2, its gravity darkening exponent (3, and its bolometric albedo A. 
Following standard practice, we set (3 = 0.25 and A = 1, which are values appropriate for 
stars with a radiative envelope. ELC uses specific intensities derived from model atmosphere 
calculations, so no parameterized limb darkening law is needed. In modeling the LMC X-1 
data, we use model atmospheres computed for the LMC metalicity from the OSTAR2002 
grid (Lanz & Hubeny 2003), supplemented by models from the BSTAR2006 grid at lower 
temperatures (Lanz & Hubeny 2007). 

ELC can include light from a fiared accretion disk, and geometrical effects due to the 
disk (e.g. eclipses of the secondary star). If the disk contributes a significant amount of 
the light from the binary the observed amplitude of the ellipsoidal light curve will decrease. 
Thus it is critical to understand how luminous the accretion disk is relative to the star. We 
show in the appendix that the accretion disk contributes at most ~ 10% of the light in the 
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optical and most likely contributes less than a few percent of the light. We therefore omit 
the disk in the models. We discuss in §4.4.21 how our results depend on the assumption of 
no disk light. 

The effects of X-ray heating are computed using a scheme adapted from Wilson (1990). 
The source of the X-rays is assumed to be a thin disk with a radius that is very small 
compared to the orbital separation. The extent to which X-ray heating changes the light 
(and velocity) curves depends primarily on how large the X-ray luminosity is compared to the 
bolometric luminosity of the secondary star. From our ongoing work on modeling 53 RXTE 
PCA X-ray spectra of LMC X-1, we find that the isotropic bolometric X-ray luminosity of 
LMC X-1 is quite steady over a time span of several years and is ~ 2.3 x lO^^erg s~^ . For 
comparison, the bolometric luminosity of the star is L ^ 1 x 10'^^ erg s~^. Because the star 
intercepts only a modest fraction of the total X-ray flux, we thus expect that the effects 
of X-ray heating on the light and velocity curves will be minor. We show in §4.4.21 that 
our main results are insensitive to reasonable changes in Lx-, and that systematic errors in 
the inclination and the derived masses due to improper treatment of X-ray heating are very 
small. In the following analysis we adopt = 2.3 x 10^^ erg s~^. 

To summarize, we have seven free parameters: i, K2, M2, P, Tq, T2, and f2- The data 
we model include 58 5-band and 57 V-band measurements from SMARTS, and 37 radial 
velocity measurements. The J-band light curve from SMARTS proved to be too noisy to 
model. The uncertainties on the individual observations were scaled to give k, N for 
each set separately. The mean of the error bars after the scaling are 0.012 mag for 5, 0.015 
mag for V", and 2.62 km s~^ for the radial velocities. Finally, we have three additional 
observational constraints that we include in the model "fitness" (see Orosz et al. 2002): the 
projected rotational velocity of the star (Vrotsin^ = 129.9 ± 2.2 km s~^), the radius of the 
secondary star {R2 = 17.0 ± 0.8 i?©), and the fact that the X-ray source is not eclipsed. 

Parameter space was thoroughly explored by running ELC's genetic optimizer several 
times with the order of the parameters changed (this gives different initial populations while 
leaving the volume of the explored parameter space unchanged). The value of T2 was confined 
to the 2a range 32, 200 < T2 < 34, 200 K as determined from the spectra, while the other 
parameters were given generous ranges. The final values of the parameters were refined using 
the grid search optimizer. The uncertainties on the fitted and derived parameters were found 
using the technique described in Orosz et al. (2002). Table [3] gives the fitted parameters, 
some derived parameters, and their uncertainties for our adopted circular orbit model and for 
an eccentric model discussed in more detail below. Figure US] shows curves of vs. various 
fitting and derived parameters of interest for the circular orbit model, from which one may 
read off the best-fitting value of a given parameter and the la, 2a, and 3a confidence ranges. 
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Each parameter is well constrained by the ellipsoidal model, with the exception of the mean 
temperature of the secondary star T2. The fact that T2 is not constrained by the ellipsoidal 
model is not surprising in this case (e.g. a single optical star in a non-eclipsing system) and 
is not a cause for concern since the optical spectra constrain the temperature quite nicely. 
We discuss in §4.4.21 how our results depend on the adopted value of T2 for a wide range of 
values. Figure [2] shows the folded light curves with the corresponding model curves, while 
the folded radial velocity curve and the corresponding model curves are shown in Figure [61 

4-4-2- Possible Systematic Errors in the Light Curve Models 

Our adopted parameters depend, among other things, an assumption that the orbit is 
circular and the star rotates synchronously with the orbit, and on two measured quantities 
that are not directly constrained (or only weakly constrained) by the light or velocity curves: 
the mean temperature of the secondary star T2 and its radius i?2- The radius of the star, 
when used as an observed quantity, helps constrain the scale of the binary and hence the 
component masses. The temperature is needed to find the radius independently. We did 
numerical experiments to see how our results change when nonsynchronous rotation and 
eccentric orbits are allowed, when the input temperature of the star is changed, and when 
the radius (used as an extra constraint) is changed. We also explored how our results depend 
on our assumption that there is no disk light contamination in the optical and how our results 
change with different assumed values of the X-ray luminosity. But first, we briefly discuss 
the slight difference between the i^-velocity derived from simple sinusoid fits to the radial 
velocities and the iiT-velocity found from the joint radial velocity and light curve fits. 

Systematic Errors in K2. The estimated i^- velocity for LMC X-1 depends slightly on 
which lines are used and how the spectra are normalized ( §3.2l and Fig. Hj). A four-parameter 
sinusoid fit to the thirty seven radial velocities measured using all of the suitable lines yields 
an orbital period of P = 3.909139 ± 0.000054 d, and K2 = 69.79 ± 0.65 km s-\ with 
= 36.50. This i^-velocity is slightly smaller than the i^-velocity found from the detailed 
modeling of the light and velocity curves {K2 = 71.61 ± 0.67 km s~^). [For comparison 
Hutchings et al. (1983) derived K2 = 65 ± 7 km s^^ for a period of P = 3.909 d.] There 
are two reasons for the difference. The main reason is that the orbital phases of the radial 
velocities are tightly constrained when they are modeled together with the light curves. In 
this case, the time of the inferior conjunction of the secondary star found from a simple 
sinusoid fit to the velocities is slightly different from the same time found from the combined 
light and velocity curve modeling (see below). The second reason for the slight difference in 
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the A'-velocities is that ELC computes corrections to the radial velocities due to the distorted 
nature of the star, so the model curve is not precisely a sinusoid, even for a circular orbit. 

Fitting the radial velocities along with the light curves provides a physical scale for the 
model. If one knew the i^-velocity ahead of time, then one in principle could fit the light 
curves by themselves. However, since the i^-velocity is not known exactly, it is better to 
use an additional constraint for K2 rather than a fixed input value. We ran fits to the 
light curves by themselves using several values of K2 as an extra constraint, assuming an 
uncertainty of 0.65 km s~^. The period was fixed at the value found from the four-parameter 
sinusoid fit, leaving six free parameters. The results are shown in Table |H Not surprisingly, 
the masses scale as K2- Since the computation of the light curves is not completely inde- 
pendent of the scale of the system (the model atmosphere specific intensities are tabulated 
in terms of the temperature and gravity in cgs units), there is a little bit of variation in 
the derived values of Tq, and . The most important result contained in Table H] is the 
fact that the derived values of the time of inferior conjunction of the secondary determined 
from the light curves is ~ 2cr different than that derived from the four-parameter sinusoid 
fit to the radial velocities. Since the photometrically determined time of conjunction should 
agree with the spectroscopically determined one, it is clearly better to fit the radial veloc- 
ities together with the light curves. The i^'-velocity derived from the simultaneous fit will 
be closer to the actual value, but at the expense of a larger value contribution from the 
radial velocities. 

Possible Nonzero Eccentricity and Nonsynchronous Rotation. Although the timescales 
for a close binary to evolve to a circular orbit and synchronous rotation are relatively short, 
we should nevertheless entertain the possibility that the orbit of LMC X-1 could be slightly 
eccentric and/or the rotation could be nonsynchronous because the age of the massive sec- 
ondary star in LMC X-1 might be less than the timescales for synchronization or circular- 
ization. Also, some examples of high mass X-ray binaries with eccentric orbits are known 
and include Vela X-1 (Rappaport, Joss, & McClintock 1976; Bildsten et al. 1997) and M33 
X-7 (Orosz et al. 2007). Having an eccentric orbit adds the eccentricity e and the argument 
of periastron uj as free parameters (see Avni 1976 and Wilson 1979 for a discussion of gen- 
eralized Roche potentials for eccentric orbits). As before, parameter space was thoroughly 
explored by using the genetic optimizer, and the results are shown in Table [31 

Avni (1976) also discussed the generalization of the Roche potential for nonsynchronous 
rotation. The figure of the star depends somewhat on how fast it is rotating, and a convenient 
parameterization of the rotation rate is fi, which is the ratio of the rotation frequency of the 
star to the orbital frequency. For a circular orbit and synchronous rotation, we have f2 = 1. 
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For an eccentric orbit, the idea of synchronous rotation is somewhat harder to define because 
the angular orbital speed of the star changes over its orbit. One can compute Q so that the 
rotation frequency of the star matches its orbital frequency at periastron (Hut 1981). We 
computed models for a wide range of Q for both circular and eccentric orbits, and the results 
are given in Table [51 

The eccentricity is formally nonzero at the ~ 3a level since the total of the fit 
decreased by about 14 with the addition of two more free parameters. However, there are 
several reasons for believing that the nonzero eccentricity is spurious: (i) The sampling of the 
radial velocities is such that a best-fit sinusoid is shifted in phase relative to the light curves. 
In an eccentric orbit, the phases of the minimum or the maximum velocity are shifted relative 
to a sine curve (see Fig. [H]), so that in the combined light-l- velocity curve modeling the total 
can be lower. Indeed, the value of the fit to the radial velocities for the eccentric 
model (Table [3]) in the combined light+velocity curve analysis is comparable to the of 
the four-parameter best-fit sinusoid fit discussed above. In other words, an eccentric orbit 
fitted to the radial velocities alone is not significantly better than a simple sine curve fit. (ii) 
The argument of periastron u is consistent with 270°. This gives one pause because tidal 
distortions of the secondary can give rise to distorted line profiles, which in turn give rise to 
systematic errors in the measured radial velocities. These velocity errors can result in fits 
with spurious eccentricities and with the argument of periastron at a quadrature phase (e.g. 
Wilson & Sofia 1976; Eaton 2008). ELC does compute corrections to the radial velocities in 
the manner of Wilson & Sofia (1976), but these corrections may not be completely accurate, 
(iii) When nonsynchronous models are considered (Table E]) we find that the eccentricity is 
correlated with the parameter Q. As the value of fl gets smaller (i.e. as the star rotates slower 
and slower), the eccentricity and the inclination get larger. In fact, values of Q smaller than 
about 0.65 are not possible since the best-fitting inclination would produce an X-ray eclipse, 
which is not observed at a high level of confidence. Hence, we would have for the best-fitting 
model in Table 5 a star rotating at ~ 65% of the synchronous value in a moderately eccentric 
orbit where the X-ray source passes just below the limb of the star as seen in the plane of the 
sky. If we ignored the fact that there is no X-ray eclipse, our best-fitting model would have 
i = 90°. (iv) The differences between the model light curves for the eccentric orbit model 
and the circular orbit model are quite subtle, (v) Detailed model computations specific to 
LMC X-1 give timescales for synchronization and circularization that are both a few orders 
of magnitude shorter than the lifetime of the secondary star for most values of the initial 
rotation rate (Valsecchi, Willems, & Kalogera, private communication 2008). 

Given these arguments which cast doubt on the reality of a nonzero eccentricity, we 
adopt the parameters derived from the circular orbit model. To be conservative with the 
uncertainties, the small differences between the parameters derived from the circular orbit 
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model and the eccentric orbit model are taken to be a measure of the systematic errors. 
These systematic errors are added in quadrature to the uncertainties found for the circular 
orbit model to produce the final adopted uncertainties (fourth column in Table [3]) • 

Regarding the models for nonsynchronous rotation, some trends are evident in Table [51 
First, as Q goes from high values to lower values, the mass of the secondary star and the 
orbital increase, and the mass of the black hole decreases (i.e. the mass ratio Q decreases). 
The reason for the mass ratio change is easy to understand. One can show that the rotational 
velocity of a star that fills its Roche lobe and is in synchronous rotation depends on the mass 
ratio (Wade & Horne 1988). The same holds true for a star that underfills its Roche lobe, 
although the function relating the velocity to the mass ratio becomes more complicated. 
When the rotation is nonsynchronous, the relation between the velocity and the mass ratio 
changes further. The reason for the change in the inclination is also easy to understand. As 
the value of Q decreases (i.e. as the star rotates more slowly), the star becomes less distorted. 
As a result, higher inclinations are needed to get roughly the same amplitudes for the light 
curves. 

Finally, as the value of Q goes down, the value of fo^' the fit goes down as well. 
However, the change in over the entire range of Q for the circular orbit is modest (Ax^ ^ 
3.5), whereas for the eccentric orbit the change is much larger (Ax^ ^ 14). We have already 
discussed why we believe that the non-zero eccentricity may be spurious. For the circular 
orbit model, we have no strong observational constraint on the value of Q. Based on other 
considerations, we give two reasons why we believe f2 = 1. First, our formal best solution 
corresponds to a near-grazing X-ray eclipse, which although possible, seems unlikely. Second, 
as discussed above, the timescales for both synchronization and circularization are much 
shorter than the age of the secondary star for most initial conditions. We conclude that Q 
is most likely to have a value very close to unity. 

Fixed Temperature. In our next experiment, we imagine that the temperature of the 
secondary is known exactly. In this case, the bolometric correction is known exactly (to 
the extent that one believes the model atmosphere computations). Once the bolometric 
correction is known, the radius of the star can be found using the apparent K magnitude, 
the distance, and the extinction Ay. We ran eight simulations where the temperature was 
fixed at values between 30,000 K and 37,000 K in steps of 1000 K. We assumed a distance 
modulus of 18.41 dz 0.10 mag and an extinction of Ay = 2.28 ± 0.06 mag. Table [6] gives 
the derived radii and luminosities for these eight cases. For each case, we fit the data using 
the genetic code, assuming the orbit is circular. Table [6] also gives the resulting values of 
the inclination i, secondary star mass M2, the black hole mass M, and the of the fit. A 
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few trends are obvious from Table [61 First, as the temperature increases, the derived radius 
is roughly constant until T = 34, 000, after which it decreases. This is due in part to the 
way the bolometric corrections change with temperature. The luminosity also increases with 
increasing temperature. Second, as the temperature increases, the of the fit generally 
goes down. On the other hand, the derived component masses do not change that much over 
the range of temperatures. The extreme values of the black hole mass are 8.79 ± 0.75 Mq 
and 11.21 ± 1.35 M^.. 



Fixed Radius. Next, we can imagine that the radius of the secondary star is known 
exactly. In this case, we do not need to know either the distance to the source or the 
extinction. Using a temperature of T2 = 33, 200 ± 500 K, we can compute the luminosity 
of the star and its uncertainty. We ran eight simulations. In each we fit the data using the 
genetic code (again using a circular orbit) with the secondary star radius fixed at one of the 
eight equally-spaced values between 15.0 i?© and 19.0 i?©. The results of this exercise (see 
Table [7]) can be used to judge what would happen to our results if our adopted distance 
to the LMC is in error, or if our adopted value of the extinction Ay were in error. For 
example, if the LMC were actually closer than what we assume, the radius that we would 
derive would be smaller (with all other things being equal). As before, a few trends are 
evident when inspecting Table [71 First, it is obvious that the luminosity of the star will go 
up as its radius goes up. As the radius increases, the scale of the binary increases leading to 
larger masses. As the star's radius goes up, it fills more of its Roche lobe, resulting in lower 
inclinations. In fact, good solutions with radii larger than ^ 19.0 M© are not possible since 
the star would then overfill its Roche lobe. Unlike the case when the temperature was fixed 
at various values, the of the fit changes very little when the radius is fixed at different 
values, except when R2 > 18 Rq. 

Variations in X-ray Heating. As noted above, ELC assumes the source of X-rays is 
a thin disk with a small diameter compared to the orbital separation. Our value of 
(Ri 2.3 X 10^^ erg s~^, which includes the geometrical factor to account for the inclination 
angle of the disk as seen from Earth) may need two correction factors when used as an 
input to the simple way in which ELC computes the effects of X-ray heating. First, if there 
is extinction local to the source (for example due to the outer edge of the accretion disk), 
some of the X-rays that would have reached the secondary star could be absorbed. If a 
flared accretion disk is present in the model, ELC can test which grid elements on the star 
are shielded from the X-ray source by the disk rim. However, it is hard to know what 
parameters to choose for the disk without more observational constraints. In any case, the 
first correction factor would be less than unity. The second correction factor comes from the 
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fact that the X-ray source may not be a perfectly flat and thin disk that is in the orbital 
plane (i.e. the disk may be warped or there may be some kind of spherical corona that emits 
some of the X-rays). When there is a perfectly thin and fiat disk, an observer on the star 
would see the X-ray emitting thin disk at an angle, and the intensity of X-rays observed 
at that point would be lower than it would be for the simple point source case owing to 
foreshortening of the disk (the foreshortening of the surface element itself is accounted for 
in the code). A warp in the disk or the presence of a quasi- spherical corona could increase 
the effective emitting area, leading to more X-rays hitting the secondary star. The second 
correction factor would be greater than unity. 

To see how much the derived parameters depend on our adopted value of L^,, we ran 
nine numerical experiments with logL^^^ fixed at values between 38.0 and 39.0 in steps of 0.25 
dex. We also set = (e.g. the case for no X-ray heating). We assume a circular orbit, 
and that the free parameters and the ranges thereof were the same as in our main work. 
Table [8] gives the results. As the value of goes up, the of the fit remains virtually 
the same. In addition, there is little variation in the inclination and in the derived masses, 
and we conclude that our overall results are rather insensitive to our treatment of the X-ray 
heating. 

Variations in the Disk Contribution. We show in the Appendix that the accretion 
disk contributes a few percent or less of the fiux in the optical bands, and at most ^ 10% 
when extreme ranges of parameters are considered. In the analysis above we assumed there 
was zero disk contribution. We ran models which included an accretion disk to see how our 
results depend on this assumption of zero disk contribution in the optical. In the ELC code 
one must specify the inner and outer disk radii, the disk opening angle, the temperature at 
the inner edge, and the exponent on the power-law temperature profile. Since LMC X-1 is 
not eclipsing, there is little to constrain the parameters of the accretion disk. Consequently 
we fixed the inner radius, the opening angle, and the exponent on the power-law temperature 
profile. The genetic code was run four times, and models where the disk fraction was within 
small threshold of a specified value were selected by means of an additional term in the total 
X^- The various runs had the disk fraction in V set at 0.05, 0.10, 0.15, and 0.20. If there is 
optical light from the accretion disk, then the 0-star by itself would be less luminous since 
the apparent magnitude of the disk+star combination is fixed. As a result its computed 
radius would be smaller for a fixed temperature. Therefore in each of the four cases the 
assumed radius of the 0-star was adjusted downward accordingly. The results are shown in 
Table [H In general, as the disk fraction increases, the inclination increases slightly and the 
masses go down slightly. For a disk fraction of 10%, the mass of the black hole decreases by 
~ lo". We therefore conclude that our results are relatively insensitive to our assumption of 



- 23 - 



no disk contribution in the optical. 



5. Discussion 

5.1. The Mass of the Black Hole 

The precision of our measurement of the black hole mass (10.91 ± 1.55 M©) represents 
a great improvement over the pioneering work of Hutchings et al. (1987) who give a mass 
"near QMq" for the black hole. Prior to this work, relatively little progress was made 
on the dynamics of LMC X-1; for example, in their recent compilation Charles & Coe 
(2006) give a range of 8 — 20 for the mass of the black hole in LMC X-1. Most of 
the binaries that contain dynamically confirmed black holes have low mass secondary stars 
(where M2 ^2.5 Mq), and the sample size of black hole binaries with high mass secondaries 
(where M2 > 10 Mq), is rather limited at the moment. There are three systems with high 
mass secondaries where dynamical studies unequivocably show that the compact object must 
be a black hole: LMC X-1, M33 X-7, and Cyg X-1. The component masses of M33 X-7 are 
well constrained (M = 15.65 ± 1.45 M© and M2 = 70.0 ± 6.9 Mq, Orosz et al. 2007), whereas 
there is some disagreement on the parameters of Cyg X-1 (Herrero et al. 1995 give masses 
of ~ 10 Mq and 18 Mq for the black hole and secondary star, respectively while Charles 
& Coe give > 4.8 for the black hole mass). There are two other eclipsing high mass 
X-ray binaries to mention in this context. 4U 1700-377 has a compact object with a mass 
of around 2.4 Mq, but this object may be a massive neutron star (Clark et al. 2002). IC 

10 X-1 has a large optical mass function that was measured using the He II A4686 emission 
line (7.64 ± 1.26 Mq, Silverman & Filippenko 2008; Prestwich et al. 2007). Assuming the He 

11 line properly tracks the motion of the Wold-Rayet (WR) secondary star, and assuming a 
mass for the WR star based on evolutionary models, the compact object in IC 10 X-1 has a 
mass of ~ 30 Mq or more. Fortunately, the future for studies of black hole binaries with high 
mass secondaries appears promising. Deep X-ray surveys of nearby galaxies are finding more 
and more high mass X-ray binaries that can be followed up in ground-based studies with the 
new generation of large telescopes and high-performance focal-plane instrumentation. 



5.2. Evolutionary Status of the Secondary Star 

Using our measured temperature and luminosity of the secondary star star (33, 200 ± 
500 K, and logL/Lg = 5.50 ± 0.05) we can place the star on a temperature-luminosity 
diagram (Figure [T^ . The position of the star in the diagram can be compared with theoret- 
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ical evolutionary tracks for single stars (Meynet et al. 1994) and with theoretical isochrones 
(Lejeune & Schaerer et al. 2001). Although we recognize that the evolutionary models for 
massive stars are uncertain owing to assumptions regarding stellar winds and rotation, for 
the purposes of this discussion we will take the evolutionary tracks and the isochrones at 
face value. The position of the LMC X-1 secondary is very close to the 5 Myr isochrone. 
Meynet et al. (1994) give evolutionary tracks for ZAMS masses of 20, 25, 40, and 60 M©. 
To get a crude idea of where the evolutionary track of a star with an initial mass of 35 Mq 
would be, we plotted the initial luminosity vs. the initial mass and the initial temperature vs. 
the initial mass for the tabulated models and used quadratic interpolation to determine the 
initial temperature and luminosity of a star with a ZAMS mass of 35 Mq. The evolutionary 
track for the star with a ZAMS mass of 40 Mq was then shifted down to the interpolated 
position of the 35 Mq star and plotted. The interpolated track with a ZAMS mass of 35 Mq 
passes very near the position occupied by the LMC X-1 secondary. The star with a ZAMS 
mass of 40 Mq has lost about 2.6 Mq by the time its temperature has fallen to ~ 33, 000 K. 
If a star with a ZAMS mass of 35 Mq behaves in a similar manner, then its mass would be 
^ 32.4 M0 when its temperature has dropped to 33,000 K, which is consistent with what we 
measure (M2 = 31.79 ± 3.48 Mq, Table E]). 

As one can see from the increasing space between the isochrones, as the age increases 
these massive stars move quickly across the temperature-luminosity diagram. Since the 
secondary star in LMC X-1 nearly fills its Roche lobe (the fraction filled is about 90%), it 
will soon enter into an interesting stage of binary star evolution. In round numbers, the 
star's current radius is 11 Rq and the radius of its Roche lobe is 19 Rq. From interpolation 
of the evolutionary track, the star with a ZAMS mass of 40 has a radius of 17 Rq at 
an age of 4.05 Myr and a radius of 19 Rq at an age of 4.22 Myr. If the evolution of a star 
with a ZAMS mass of 35 Mq proceeds at a similar pace, the secondary star in LMC X-1 will 
encounter its Roche lobe in only a few hundred thousand years. Since the secondary star is 
substantially more massive than the black hole, the resulting mass transfer will be rapid and 
may be unstable (Podsiadlowski, Rappaport, & Han 2003). If the mass transfer is stable 
and the binary survives the initial phase of rapid mass transfer, the present-day secondary 
may become detached again in a binary with a longer orbital period. On the other hand, if 
the rapid mass transfer is unstable, the orbit will rapidly shrink and the binary will enter 
a "common envelope" (CE) phase. The outcome of the CE phase depends on how tightly 
bound the envelope of the star is. If the envelope is tightly bound, then the energy liberated 
as the two stars move towards each other will not be enough to expel the envelope, and the 
two stars will merge. On the other hand, if the envelope is loosely bound, then the envelope 
of the star can be expelled before the cores merge, leaving behind a tight binary consisting of 
the present-day black hole and the core of the present-day secondary star. Although detailed 
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computations are required to assess these scenarios, it seems that the most hkely outcome of 
the CE phase in LMC X-1 would be a merger since the envelopes of massive stars are tightly 
bound, and the density gradient in a main sequence star is not nearly as steep as it is in a 
giant (e.g. Podsiadlowski et al. 2003). What the merger product would look like is an open 
question. 



5.3. The O-star Wind and the ASM X-ray Light Curves 



The X-ray minimum/maximum in the ASM X-ray light curves occurs at inferior/superior 
conjunction of the secondary star (Fig. [H]), as expected if the modulation is caused by scat- 
tering or absorption in a stellar wind. A simple cosine wave gives a reasonable fit to the 
light curves; the respective full amplitudes for the A, B and C channels are 0.072 ± 0.010, 
0.077 ± 0.11, and 0.038 ± 0.029 (Fig. M and Table [I]). Thus, the modulation is largely in- 
dependent of energy, which indicates that electron scattering is the main source of opacity. 
We now show that a simple model of Thomson scattering in a spherically symmetric stellar 
wind is not only reasonably consistent with the observed X-ray modulation but also yields 
plausible values of the wind mass flux, accretion rate, and accretion efficiency. For this sec- 
tion we use the values of component masses and other system parameters in the right-hand 
column of Table [31 

We assume that the O star has a radiatively driven wind like those from other hot 
stars and that the wind velocity has a radial dependence that follows a "beta law", i.e., 
v{r) = Vr^(l — h/rY where b is approximately equal to the O-star photospheric radius R2 
(e.g.. iDessart fc Owockill2005l : iKudritzki fc Pulsll2000l . and references therein). The terminal 



velocity v^o can be es timated from the escap e velocity at the stellar surface and the photo- 
spheric temperature (IKudritzki fc Pulsll2000l ). The escape velocity fesc can be estimated in 
turn from the mass and radius results in Table [31 We obtain v^sc ~ 590 ± 40 km s~^. In 
the calculations below, we take the values of Voo/vqsc and of (3 to be like those in the wind 
model presented for M33 X-7 by Orosz et al. (2007), and so we adopt the values foo = 1400 
km s~^ ~ 2.4fesc and (3 = 1. This valu e of Vnc^/vesr. is appropri ate for an O-star with pho- 
tospheric temperature T2 ~ 33,000 K (IKudritzki fc Pulsll2000l ). The wind velocity at the 
radius corresponding to the location of the black hole (r = a = 36.5 Rq) is then v{a) ~ 740 
km s~^. Adding in quadrature the transverse velocity of the black hole through the wind, we 
find that the gas fiows by the black hole at a velocity V ~ 880 km s~^. Given the substantial 
uncertainties in both the form of the velocity law and the values of the parameters, these 
estimates of v{a) and V must each be uncertain by as much as 50%. 



In any spherically-symmetric wind model, the scattering optical depth along the line 
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of sight to the X-ray source at superior conjunction only differs from that at inferior con- 
junction by the optical depth along an additional path of length d = 2a sin i = 43.3i?0 at 
superior conjunction. Let be the electron column density along this part of the superior 
conjunction line of sight. Then there must be a scattering optical depth along this part of 
the line of sight of NcCTt ~ 0.07, where ctt is the Thomson cross section. This, in turn, 
implies iVe ~ 1.1 x 10^^ cm~^. The mean electron number density along the path of length 
d is then He ~ 3.5 x 10^° e~ cm~'^. In our adopted wind model, this mean density is ~ 33% 
higher than the density at the black hole orbital radius (and it is also ~ 33% lower than at 
the point of highest density along this section of the line of sight). The gas density at the 
black hole orbital radius is then p{a) ~ 5.3 x 10~^^ g cm~^. The corresponding mass loss 
rate in the 0-star wind is Mw = 4na'^p{a)v{a) ~ 5 x 10~^ MQ/yr~-^. 

If the wind of the O star in LMC X-1 is only partially ionized like the winds of ordinary 
0-stars, then photoelectric absorption should be manifest via a larger degree of modulation 
in the lower energy band ASM light curves. Is it a reasonable assumption in the case of 
LMC X-1 that iron and the other metals in the wind have been fully stripped of their K 
shell electrons? Some support for this idea is provided by the remarkable photoionized He 
III nebula that surrounds LMC X-1 (PakuU & Angebault 1986). More directly, at every 
point on the path described above the ionization parameter (Tarter et al. 1969) is large, 
L/nB? > 1000. Furthermore, the spectrum of LMC X-1 is soft (disk blackbody temperature 
kT ^ 0.9 keV), which makes it an effective ionizing agent. Under these conditions, studies 
of similar high-mass X-ray binaries support our conclusion that the wind is close to fully 
ionized in the region in question (McCray et al. 1984; Vrtilek et al. 2008). 

A crude estimate of the accretion rate is given by the rate at which the matter in the 
wind enters a cylinder of radius ~ 2GM/V'^ centered on the black hole, where G is the 
gravitational constant, M is the black- hole ma ss, and V is the (above est imated) velocity 



of the wind relative to the black hole (see, e.g., IShapiro fc Teukolskylll983l ). The accretion 
rate is then given by Mb_h ~ ^'"capP('^)^- Using the value of M in Tableland the values 
of p(a) and V estimated above, we obtain Tc ~ 5 Rq and Mb-h ~ 3 x 10"*^ Mq jt~^ with 
uncertainties in each value of roughly a factor of two. 

As noted above, we find that the isotropic bolometric X-ray luminosity of LMC X-1 is 
quite steady over a time span of several years and is ~ 2.3 x 10^^ erg s~^ . With M = 10.91 
Mq and the above estimated accretion rate, the implied accretion efficiency is ?7 ~ 0.1. This 
is comparable to the 0.06 efficiency of a Schwarzschild black hole or the canonical value of 
0.1 that is commonly used, but, considering that this efficiency estimate is uncertain by at 
least a factor of two, this may be fortuitous. 



6. Summary 



We present a new dynamical model of the high mass X-ray binary LMC X-1 based 
on high quality echelle spectra, extensive optical light curves, and infrared magnitudes and 
colors. From our optical data we find an orbital period of P = 3.90917 ± 0.00005 days. 
We present a refined analysis of the All Sky Monitor data from RXTE and find a period of 
P = 3.9094 ± 0.0008 days, which is consistent with the optical period. A simple model of 
Thomson scattering in the stellar wind accounts for the 7% modulation seen in the X-ray light 
curves. We find that the V^-band extinction to the source is much higher than previously 
assumed. The V — K color of the star (1.17 ± 0.05) implies Ay = 2.28 ± 0.06, whereas 
several determinations of the column density from X-ray observations give Ay > 2.57 ±0.12. 
The color excess of E{B — V) = 0.44 together with Ay = 2.27 gives Ry = 5.18, which is 
much higher than the nominal mean value in the LMC of Ry = 3.1. For the secondary 
star we measure a radius of R2 = 17.0 ± 0.8 i?0 and a projected rotational velocity of 
Vj-otsini = 129.9 ± 2.2 km s^^. Using these measured properties of the companion star to 
constrain the dynamical model of the light and velocity curves, we find an inclination of 
i = 36.38 ± 1.92°, a secondary star mass of M2 = 31.79 ± 3.48 M©, and a black hole mass 
of 10.91 ± 1.41 Mq. The present location of the secondary star in a temperature-luminosity 
diagram is consistent with that of a star with an initial mass of 35 Mq that is 5 Myr past 
the zero-age main sequence. The star nearly fills its Roche lobe (~ 90%) and, owing to the 
rapid change in radius with time in its present evolutionary state, it will encounter its Roche 
lobe and begin rapid and possibly unstable mass transfer in only a few hundred thousand 
years. 
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A. Accretion Disk Contamination in the Optical and Infrared Bands 

The optical and infrared bands of interest to us for the elhpsoidal modeUng (B, V, and J) 
are in the Rayleigh- Jeans part of the spectrum for the emission from both the secondary star 
and the disk (see below). Therefore, the flux we receive in any particular band is proportional 
to the product of the projected area of the radiating surface A and the temperature T. 

The radius of the secondary star is 17.0 Rq and its photospheric temperature is 33, 200 K. 
Therefore, for the radiation from the star, we have 

(Ar)sec = 1.5 X 10^^ cm^ K. (Al) 

In the case of the disk, we note that the circularization radius of the captured wind material 
is likely to be quite small. Nevertheless, the steady state disk will probably have a large 
outer radius since the gas needs to get rid of its angular momentum. An extreme possibility 
is that the angular momentum is removed solely by tidal torques from the companion star 
acting on the outer regions of the disk. From the discussion in Frank, King & Raine (2002; 
eq. 5.122), the outer radius of the disk is then approximately equal to 0.9 times the Roche 
lobe radius. Using a standard approximation for the Roche lobe radius (Frank et al. 2002), 
we find for the outer radius 

/ Mrh V^^ 

i?out = 0.9 X 0.462 — a = 6.6 X 10^^ cm. (A2) 

We further note that most of the disk will be dominated, not by local viscous dissipa- 
tion, but by reprocessing of radiation emitted by the inner region of the disk. Assuming 
that 10,^10% of the bolometric luminosity of the disk is reprocessed (the exact fraction is 
uncertain), the luminosity of the reprocessed radiation from the disk is 2.3 x lO^^^io ergs^^. 
For a typical reprocessing geometry, the effective temperature of the reprocessed radiation 
will vary as (as compared to the steeper R^^^^ variation one expects for viscously 

generated flux). Let us write the disk temperature as 

T{R) = Tio (i?/10^° cm)-i/^ (A3) 

and let us assume that the reprocessing region of the disk extends from i?in ~ lOORs — 
3.1 X 10^ cm to i?out (the precise value of i?in is not important for what follows). Then we 
have the condition (in cgs units) 

/•6.6X10" 

2.3 X lO^^Cio = / iirRaT^oW^^^^y^dR. (A4) 

^3.1x108 

Solving, we obtain 

rio = 8.1xl0^ei'o^'K. (A5) 
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The disk temperature thus varies from 4.6 x lO^^jg^ K at to 1.0 x lO^^Jg^ K at -Rout- Even 
the outermost region of the disk is hot enough that we may treat its optical/IR emission in 
the Rayleigh- Jeans hmit. 

We can now estimate the value of (AT) for the disk emission: 

;>6.6xl0" 

(AT)disk = / 27ri?(cosi)Tio {R/lO^^y^^^dR = 1.5 x lO'^^Cw^ cm^ K. (A6) 

Comparing this estimate with equation (lAip . we see that the disk will contribute about 10% 
of the flux in the optical and infrared. Note that the uncertain parameter appears only 
with a 1/4 power, so we are not sensitive to its value. However, the result is very sensitive 
to the assumed outer radius of the disk. If the disk is smaller than we have assumed, the 
disk flux will be reduced substantially, roughly as -Rout- 

For a wind- fed system like LMC X-1, the disk will probably get rid of much of its 
angular momentum by interacting directly with the incoming material. The efficiency of 
this process is very difficult to estimate from first principles, so we will use the eclipsing 
system M33 X-7 as a guide to what the outer radius of the accretion disk might be. Orosz 
et al. (2007) determined an outer radius of 0.45 ± 0.03 times the Roche radius for M33 X-7. 
If the accretion disk in LMC X-1 fills 45% of its Roche lobe then the disk contribution will 
be (v4T)disk = 5 X 10^^,^}o^, which is about 3.5% of the optical/infrared fiux. Even this may 
be an overestimate. It is quite conceivable that the disk is much smaller than 45% of the 
Roche lobe radius, which would make the disk contamination completely insignificant. 

We have ignored reprocessing of radiation from the secondary star, but one can show 
that it is again not important. 
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Table 1. Fitting Parameters for ASM Light Curves'^ 



Channel 




ai 


02 






A 


0.7291 ± 0.0025 


-0.0263 ±0.0036 


-0.0031 ± 0.0035 


1.942 


0.072 ± 0.010 


B 


0.5275 ± 0.0020 


-0.0204 ±0.0028 


0.0037 ±0.0028 


1.151 


0.077 ±0.011 


C 


0.2602 ± 0.0026 


-0.0050 ±0.0038 


0.0005 ± 0.0037 


0.681 


0.038 ± 0.029 


SUM 


1.5201 ± 0.0042 


-0.0491 ± 0.0060 


-0.0004 ± 0.0059 


1.236 


0.065 ± 0.008 



= 00 + 01 cos(27r0) + 02 sin(27r0) 
''fractional amplitude A= (2|ai|/ao) 



Table 2. LMC X-1 Column Density Measurements 



Model Number 


Model 


Mission 


Nh (102icm-2) 


Av (mag) 


Reference 


1 


DBB 


Chandra HETG 


4.6 ±0.2 


2.57 ±0.12 


Cui et al. 2002 


2 


COMPTT 


Chandra HETG 


5.9 ±0.3 


3.30 ±0.18 


Cui et al. 2002 


3 


DBB 


BeppoSAX 


8.1 ±0.5 


4.53 ±0.29 


Haardt et al. 2001 


4 


GRdisk 


ASCA 


5.3 ±0.2 


2.96 ±0.12 


Gierlinski et al. 2001 


5 


DBB 


ASCA 


6.3 ±0.1 


3.52 ± 0.08 


Nowak et al. 2001 


6 


DBB 


BBXRT 


5.8 ±0.9 


3.24 ±0.51 


Schlegel et al. 1994 


7 


X-ray halo 


Chandra 


6.5 ±0.1 


3.65 ± 0.10 


Xiang et al. 2005 
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Table 3. LMC X-1 Adopted Parameters 



parameter 


value 


value 


Adopted 




(circular orbit) 


(eccentric orbit) 


value*^ 


P (days) 
Tq (HJD 2,453,300+)'' 


3.90914 ± 0.00005 


3.90917 ± 0.00005 


3.90917 ± 0.00005 


91.3436 ± 0.0078 


91.3419 ± 0.0087 


91.3436 ± 0.0080 


A2 (km s ) 


71. ol ± U.o7 


70.74 ± 0.75 


71.61 ± I.IU 


i (deg) 


Qo _i_ 1 no 
00.00 ± i.yz 


/ .UU ± i. 10 


OP. OQ _i_ no 
00.00 ± Z.KjZ 




0.886 ±0.036 


0.894 ±0.040 


0.886 ± 0.037 


M2 (Mq) 


31.79 ±3.48 


30.62 ±3.17 


31.79 ±3.67 


e 




0.0256 ± 0.0066 




a; (deg) 




260.5 ± 16.8 




To (HJD 2,453,300+)^ 




91.3072 ± 0.2073 




i?2 (-Rq) 


16.89 ±0.87 


16.41 ± 0.71 


17.00 ± 0.80*^ 


logg (cgs) 


3.485 ±0.014 


3.497 ±0.011 


3.485 ±0.018 


a (i?©) 


36.49 ± 1.42 


35.97 ± 1.28 


36.49 ± 1.51 


M (Mo) 


10.91 ± 1.41 


10.30 ± 1.18 


10.91 ± 1.54 


(5, \/ bands) 


59.18, 60.44 


59.00, 59.70 




(radial velocities) 


49.63 


36.20 




(total) 


169.27 


155.39 





Note. — The effective temperature of the secondary has been constrained to the 
range 32, 200 < T^s < 34, 200 K. 

'^The uncertainties include systematic errors 

^Time of inferior conjunction of secondary 

^Time of periastron passage of secondary 

'^Determined from the temperature, K magnitude, distance, and extinction 



-36- 



Table 4. LMC X-1 Parameters as a Function of K2 



Assumed K2 


i 




M2 


M 


Xmin 


(km s~-^) 


(deg) 




(Mo) 


(Me) 




68.50 ±0.65 


36.20 ±2.02 


0.0111 ±0.0055 


31.46 ±3.51 


10.34 ± 1.31 


115.68 


69.50 ±0.65 


36.17 ±2.02 


0.0111 ±0.0057 


31.72 ±3.53 


10.58 ± 1.33 


115.66 


69.80 ±0.65 


36.17 ±2.03 


0.0111 ±0.0055 


31.79 ±3.48 


10.65 ± 1.33 


115.65 


70.50 ±0.65 


36.16 ±2.03 


0.0111 ±0.0055 


32.00 ±3.49 


10.82 ± 1.36 


115.63 


71.50 ±0.65 


36.23 ±2.01 


0.0112 ±0.0055 


32.10 ±3.51 


11.00 ± 1.38 


115.60 


72.50 ±0.65 


36.19 ±2.02 


0.0112 ±0.0055 


32.37 ±3.51 


11.26 ± 1.41 


115.57 



^A(f) is the phase shift of the photometric Tq relative to the spectroscopically deter- 
mined value. 



Note. — The assumed distance is 18.41 ± 0.10 mag, the assumed temperature is in 
the range 32, 200 < T^s < 34, 200 K, and the assumed extinction is Ay = 2.28 ± 0.06 
mag. 
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Table 5. LMC X-1 Parameters for Nonsyncfironous Rotation 



11 


i 


e 


UJ 


IVI2 


M 


2 

A,min 








(deg) 








1.15 


32.0 ± 1.0 


(fixed) 




30.7 ±2.3 


12.5 ± 1.0 


171.53 


1.10 


33.0 ± 1.4 


(fixed) 




31.4 ±2.4 


12.1 ± 1.2 


170.45 


1.05 


34.9 ± 1.5 


(fixed) 




31.0 ± 1.5 


11.3 ± 1.4 


169.91 


1.00 


36.4 ± 1.9 


(fixed 




31.2 ±3.5 


10.9 ± 1.4 


169.38 


0.95 


38.5 ±2.2 


(fixed) 




32.2 ±3.1 


10.3 ± 1.0 


169.11 


0.90 


42.8 ±2.5 


fixed 




31.6 ±2.9 


9.3 ± 1.1 


168.73 


0.85 


44.5 ±2.6 


(fixed) 




32.3 ± 2.9 


8.4 ±0.9 


168.51 


0.80 


48.1 ±2.6 


(fixed) 




33.0 ± 2.2 


8.4 ± 1.1 


168.18 


0.75 


52.1 ±2.5 


(fixed) 




33.8 ±3.1 


7.9 ±0.9 


167.84 


0.70 


56.3 ±2.7 


(fixed) 




35.6 ±2.9 


7.7 ±0.7 


168.17 


0.65 


60.9 ± 1.8 


(fixed) 




37.4 ± 2.0 


7.5 ± 0.4 


168.06 


1.15 


33.8 ± 1.5 


0.0215 ± 0.0038 


267.0 ± 17.2 


28.2 ±2.9 


11.0 ±1.2 


162.45 


1.10 


34.6 ± 1.7 


0.0210 ±0.0063 


265.5 ± 16.4 


29.4 ± 3.2 


11.0 ±1.3 


160.46 


1.05 


36.0 ± 1.8 


0.0245 ±0.0036 


261.3 ± 16.8 


29.6 ±3.3 


10.7 ± 1.4 


157.03 


1.00 


37.4 ±2.1 


0.0269 ±0.0038 


266.3 ± 16.9 


30.3 ±3.4 


10.3 ± 1.4 


154.93 


0.95 


39.4 ±2.0 


0.0287 ±0.0048 


266.3 ± 17.2 


31.4 ±3.2 


9.8 ± 1.1 


153.19 


0.90 


42.1 ± 1.7 


0.0293 ±0.0042 


264.5 ± 16.6 


32.2 ±2.4 


9.2 ±0.9 


152.35 


0.85 


44.2 ±2.2 


0.0311 ±0.0046 


265.2 ± 17.0 


33.5 ±2.8 


9.0 ± 1.2 


150.76 


0.80 


47.6 ±2.6 


0.0336 ± 0.0045 


266.2 ± 17.5 


34.5 ± 3.6 


8.6 ± 1.1 


149.94 


0.75 


50.3 ±3.5 


0.0360 ±0.0038 


268.3 ± 17.1 


36.4 ± 3.5 


8.4 ± 1.0 


149.48 


0.70 


56.8 ±4.0 


0.0372 ± 0.0035 


268.2 ± 16.5 


36.1 ±3.9 


7.5 ± 0.9 


148.32 


0.65 


61.7 ±2.6 


0.0372 ± 0.0037 


259.9 ± 17.0 


38.0 ± 2.4 


7.4 ± 0.4 


148.46 



^Q, is the ratio of the rotation frequency of the star to the orbital frequency. 

Note. — The temperature is in the range 32, 200 < T^s < 34, 200 K, the assumed 
distance is 18.41 ± 0.10 mag, and the assumed extinction is Ay — 2.28 ± 0.06 mag. 
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Table 6. LMC X-1 Parameters as a Function of Temperature 



temperature 


derived radius 


logL 


i 


M2 


M 


Xmin 


(K) 


(Re) 




(deg) 


{Me) 


{Me) 




30,000 


17.03 ±0.83 


5.33 ±0.04 


36.29 ±2.10 


32.34 ± 3.48 


11.09 ± 1.34 


174.14 


31,000 


16.89 ±0.82 


5.38 ±0.04 


36.47 ±2.11 


31.93 ± 1.6 


10.95 ± 1.34 


174.66 


32,000 


17.07 ±0.83 


5.44 ± 0.04 


36.12 ±2.03 


32.73 ±3.56 


11.21 ± 1.35 


173.94 


33,000 


17.07 ±0.82 


5.49 ± 0.04 


36.30 ±2.02 


32.45 ± 3.73 


11.07 ± 1.27 


171.30 


34,000 


16.72 ±0.81 


5.53 ±0.04 


37.11 ±2.17 


30.67 ±3.61 


10.44 ± 1.32 


169.24 


35,000 


16.23 ±0.78 


5.55 ±0.04 


38.42 ±2.33 


28.66 ±3.21 


9.66 ± 1.25 


169.55 


36,000 


15.79 ±0.76 


5.58 ±0.04 


39.92 ±2.13 


27.12 ±3.00 


8.96 ±1.13 


170.07 


37,000 


15.42 ± 0.74 


5.61 ±0.04 


40.16 ± 1.71 


26.81 ± 1.91 


8.79 ±0.75 


165.82 



Note. — The assumed distance is 18.41±0.10 mag and the assumed extinction is Ay — 2.28±0.06 
mag. 
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Table 7. LMC X-1 Parameters as a Function of Radius 



radius 


derived luminosity 


i 








M 


Xmin 


(Re) 


(log L/Lq) 


(deg) 






(Me) 


{Me) 




15.0 


5.39 ±0.03 


41.24 ±0, 


.87 


24, 


.56 ±0.95 


8.15 ±0.25 


168.52 


15.5 


5.42 ± 0.03 


39.64 ± 0, 


.80 


26. 


.34 ±0.90 


8.86 ±0.28 


168.71 


16.0 


5.45 ± 0.03 


38.17 ±0. 


,81 


28. 


.17 ±0.97 


9.62 ±0.29 


168.99 


16.5 


5.47 ±0.03 


36.82 ± 0. 


,71 


30. 


.16 ±0.89 


10.43 ±0.28 


169.19 


17.0 


5.50 ±0.03 


35.62 ±0. 


,71 


32. 


,28 ±0.97 


11.29 ±0.31 


169.44 


17.5 


5.53 ±0.03 


34.42 ± 0. 


,70 


34. 


,50 ± 1.24 


12.22 ±0.30 


169.77 


18.0 


5.55 ±0.03 


33.37 ±0, 


,65 


36. 


,86 ±0.87 


13.20 ±0.37 


170.42 


18.5 


5.57 ±0.03 


32.39 ±0, 


,56 


39. 


,29 ±0.59 


14.22 ±0.34 


170.36 


19.0 


5.56 ±0.03 


32.48 ± 0. 


,59 


42. 


,60 ±0.56 


15.44 ±0.37 


171.27 



Note. — The assumed temperature is in the range 32, 200 < Teflf < 34, 200 K, the 
assumed distance is 18.41±0.10 mag, and the assumed extinction is Ay — 2.28±0.06 
mag. 
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Table 8. LMC X-1 Parameters as a Function of X-ray Heating 



logL^ (erg s ^) 


i 


M2 


M 


Xmin 


(isotropic equivalent) 


(deg) 


(Me) 


(Me) 




no heating 


36.02 ± 1.78 


31.60 ±3.11 


10.98 ± 1.28 


169.57 


38.00 


36.02 ± 1.82 


31.61 ±3.03 


10.99 ± 1.40 


169.57 


38.25 


36.01 ± 1.28 


31.54 ±2.27 


10.99 ±0.99 


169.39 


38.50 


36.06 ± 1.83 


31.41 ±3.28 


10.95 ± 1.40 


169.35 


38.75 


36.02 ±1.39 


31.40 ±2.50 


10.97 ±0.77 


169.48 


39.00 


35.91 ± 1.77 


31.49 ±3.17 


11.05 ± 1.16 


170.30 



Note. — The assumed temperature is in the range 32, 200 < T^s < 
34, 200 K, the assumed distance is 18.41 ±0.10 mag, and the assumed extinc- 
tion is Av = 2.28 ± 0.06 mag. 



Table 9. LMC X-1 Parameters as a Function of the Disk Fraction 



Disk fraction 


derived radius 


i 


M2 


M 


Xmin 


(y-band) 


(Re) 


(deg) 


(Me) 


(Me) 




0.05 


16.58 ±0.80 


37.31 ± 1.84 


28.88 ±2.56 


10.03 ± 1.10 


169.78 


0.10 


16.15 ±0.80 


37.99 ±2.28 


27.25 ±2.65 


9.52 ± 1.12 


170.16 


0.15 


15.69 ±0.80 


38.44 ± 1.43 


25.96 ± 1.89 


9.15 ±0.73 


171.04 


0.20 


15.22 ±0.80 


39.10 ± 1.36 


24.43 ± 1.70 


8.69 ±0.60 


175.73 



Note. — The assumed temperature is in the range 32, 200 < T^s < 34, 200 K, the 
assumed distance is 18.41 ±0.10 mag, and the assumed extinction is Ay — 2.28 ±0.06 
mag. 
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Fig. 1. — 1x1 arcminute finding charts for tlie field surrounding LMC X-1 obtained with 
the Magellan-Baade telescope. The left panel is a y-band image obtained with the IMACS 
imaging spectrograph, sampled at 0'.'22 pixel"^ under seeing of O'.'S. The right panel is the 
same field observed in the i^-band using the PANIC camera with 0'.'125 pixel"^ sampling 
and a seeing of 0'.'55. The correct counterpart together with the nearby B5 supergiant R148 
are marked. 
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Fig. 2. — The phased hght curves of LMC X-1. Phase zero corresponds to the time of the 
inferior conjunction of the secondary star. Shown from the top are the B, V, J hght curves 
with the best-fitting eUipsoidal models assuming a circular orbit (solid lines) and assuming 
an eccentric orbit (dashed lines). Owing to the large scatter, the J band data were not used 
in the modeling. 
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Fig. 3. — The average spectrum of LMC X-1 in the restframe of the secondary star (black 
dots) obtained with the MIKE echelle. Most of the strong features are labeled, and inter- 
stellar lines are marked with 'IS'. The best-fitting model (red line) has T^g = 33,225 K, 
log (? = 3.56, and sini = 128.0 km s~^. The noise near 5060 A is caused by the transition 
from the blue arm to the red arm. 
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Fig. 4. — Summary of fitted values of the velocity semiamplitude K (top) and the time of 
maximum velocity Tmax (bottom) for different combinations of bandpasses, template spectra 
and modes of filtering. The lines or line group are indicated at the bottom of the figure 
(see text). The open and filled circles are for template HD 93843 and Legendre filtering. 
Triangles: HD 101205, Legendre filtering. Squares: Synthetic template, Legendre filtering. 
Crosses: Synthetic template, Fourier filtering. The filled circles and dashed hnes correspond 
to the values of K and Tmax given in the text, which were computed for the mean velocities. 
The uncertainties shown are purely statistical, correspond to = 1, and are at the \a 
level of confidence. No results are shown for the He I A4471 line for HD 101205 because the 
template spectrum was corrupted in this band. 
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Fig. 5.— The equivalent widths of HIO (top), He I A4471 (center), and He H A4200 (bottom) 
as a function of orbital phase, where phase zero corresponds to the time of the inferior 
conjunction of the secondary star. 
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Fig. 6. — Top: The phased velocity curve of LMC X-1. Phase zero corresponds to the time 
of the inferior conjunction of the secondary star. The MIKE radial velocities are shown with 
the filled circles and the MagE radial velocities are shown with the filled triangles. The model 
curve for the best-fitting orbital model assuming a circular orbit is shown with the solid line, 
and the model curve for the best-fitting model assuming an eccentric orbit is shown with the 
dashed lines. Center: The residuals with respect to the circular orbit model. Bottom: The 
residuals with respect to the eccentric orbit model. 
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Fig. 7. — Top: computed using a three-parameter sinusoid vs. the trial period for the 
Magellan radial velocities. Center: The periodogram (computed in a similar manner as 
above) for the radial velocities from Hutchings et al. (1983, 1987) combined with the Magellan 
radial velocities. Bottom: The periodogram derived from the SMARTS photometry and the 
Magellan radial velocities (see text). The best-fitting period is found to be 3. 90917± 0.00005, 
and all other possible alias periods are ruled out at high confidence. In each of the three 
panels the X-ray period given in Levine & Corbet (2006) is indicated with an open circle, 
and the refined X-ray period is denoted by the filled circle. 
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Fig. 8.— Power density spectrum (PDS) of a 1.5-12 keV ASM light curve of LMC X-1 that 
has been modified to remove variabiUty on time scales longer than ?a 30 days (sec text). The 
original FFT was oversampled by a factor of four and had approximately 62,000 frequencies 
from cycles day^^ to the Nyquist frequency of 3.3 cycles day~^. Top: The low frequency 
part of a rebinned PDS in which the number of frequency bins was reduced by a factor of 
10 by using the maximum power in each contiguous set of 10 frequency bins in the original 
PDS as the power of the corresponding bin in the rebinned PDS. Bottom: A portion of the 
original PDS. The horizontal bar above the peak shows our estimate of the possible values of 
the centroid frequency, i.e., 0.25580 ± 0.00005 d~^. In both panels, the power is normalized 
relative to the PDS-wide average. 
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Phase (P = 3.9094 days) 

Fig. 9.— The RXTE ASM light curves of LMC X-1 folded using a period of 3.9094 days 
and an epoch of phase zero of MJD 53390.75174. The four panels from bottom to top show 
the data from the A (1.5-3 keV), B (3-5 keV), and C (5-12 keV) bands, and the sum of 
the three bands (nominally comprising photon energies 1.5-12 keV). The error bars indicate 
±lcr statistical uncertainties. The smooth curves are the best-fit functions of the form 
/(0) = ao + ai cos(2'7r0) + 02 sin(27r0) (see Table [1] for the best-fitting parameters). 
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Fig. 10. — The average profiles of the four hehum fines used to determine tfie rotational 
velocity ("histograms"), and the best-fitting model profiles for each (smooth curves). The 
line identification and the derived value of l^ot sin i in km s~^ are given in each panel. 
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Fig. 11. — The radius of the companion as a function of the extinction Ay for the \/-band 
(triangles) and the i^'-band (circles). The radius is the same for both bands when the 
extinction is ~ Ay = 2.2 mag. 
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Fig. 12. — The average lUE spectrum of LMC X-1, dcrcddcncd using Ay = 2.27 and Ry = 
5.16, is shown with a model spectrum from the OSTAR2002 grid having Tcff = 32,500 K 
and logs' = 3.5, scaled using a distance modulus of 18.41 and a stellar radius of 17.08 i?©. 
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Fig. 13. — Curves of vs. various fitting and derived parameters of interest for the circular 
orbit model. The horizontal dashed lines denote the 1, 2, and 3a confidence limits. 
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log(T,ff) (K) 

Fig. 14. — The solid lines show evolutionary tracks for an LMC metahcity on a temperature- 
luminosity diagram for stars with ZAMS masses of 60 M©, 40 Mq, 25 Mq, 20 Mq taken from 
Meynct ct al. 1994). The track for a ZAMS mass of 35 Mq was crudely interpolated. The 
dotted lines show isochrones (also for an LMC metalicity, Lejeune & Schaerer 2001) from 
left to right of 0, 1, 2, 3, 4, 5, and 6 Myr. The location of the LMC X-1 secondary is very 
close to the interpolated track for the 35 Mq star and the 5 Myr isochrone. 



